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Abstract 



The chiral phase transition at high temperature is investigated using the ef- 
fective potential in the framework of the QCD-like gauge theory with a vari- 
ational approach. We have a second order phase transition at T c = 136MeV. 
We also investigate numerically the temperature dependence of (qq), f n and 
a2(T) (coefficient of the quadratic term in the effective potential) and estimate 
the critical exponents of these quantities. 
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I. INTRODUCTION 



The properties of Quantum Chromodynamics(QCD) at finite temperature are a subject 
of general interest. At zero temperature and zero chemical potential, the chiral symmetry 
is spontaneously broken. On the other hand, we expect that the chiral symmetry will be 
restored at sufficiently high temperature and it may be realized in high energy heavy-ion 
collisions at RHIC and LHC. The restoration of the chiral symmetry is important also in 
physics of neutron stars and early universe. 

Effective theories of QCD such as a model |l| and Nambu-Jona-Lasinio model or 
lattice QCD |3j have been used in the study of the chiral phase transition and the related 
critical phenomena in QCD. An effective potential and/or Schwinger-Dyson equation for 
quark propagator have been widely used to study the spontaneous chiral symmetry breaking 
and its restoration at finite temperature and/or finite chemical potential. 

In this paper we study the chiral phase transition at finite temperature in the chiral limit. 
We use the Cornwall- Jackiw-Tomboulis(CJT) type effective potential [|J for a composite 
operator in the framework of the renormalization group improved gauge theory(the so called 
QCD-like gauge theory) [5,6,7]. We take into account a correct asymptotic behavior of the 
quark self energy and use the variational method in calculating the effective potential. The 
present approach does not break the chiral symmetry explicitly in contrast to some previous 
studies of the effective theories of QCD. Using the effective potential we study the critical 
temperature and order of the chiral phase transition. If the transition is second order, we 
can calculate the critical exponents. Rajagopal and Wilczek have studied the second order 
phase transition using three dimensional 0(4) Heisenberg model |§. However it is still an 
open question whether QCD with two massless quarks and 0(4) Heisenberg model belong 
to the same universality class or not. The critical exponents have been also studied using 
lattice QCD §. 

Furthermore, in previous studies of the effective theories which respect the 
QCD[9,10,11,12] it is not obvious whether the critical exponents agree with those of the 
mean field approximation a la Landau or not. We numerically study the critical exponents 
of the order parameter for the chiral phase transition (qq)r, the pion decay constant at finite 
temperature f n {T) and a 2 (T) which is the coefficient of the quadratic term in the effective 
potential. 

This paper is organized as follows. In Sec. 2 we briefly review the effective potential 
for a composite operator in the QCD-like gauge theory. In Sec. 3 we calculate the effective 
potential at finite temperature. In Sec. 4 we present numerical results and determine the 
temperature dependence of (qq) and /„.. We also estimate the critical exponents. Sec. 5 
contains the summary and discussion. 

II. EFFECTIVE POTENTIAL AT ZERO TEMPERATURE 

In this paper we fix the mass scale by the condition Aqcd — 1- The CJT effective 
potential for the quark propagator takes the following form in Minkowski space: 

V[G,D] = J tr (HS-\p)G(p)) - S-\p)G(p) + l) 



2 



iff dp d q 



2 J J (2ir)H (2tt)H 
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where C 2 = (A" 2 — 1)/(2N C ) is the quadratic Casimir operator for color SU(N C ) group, S(p) 
is the propagator of a massless quark, and D^ij)) is the gluon propagator in the Landau 
gauge: 

D^(p) = -(g^-^-), 

which is diagonal in the color space and "tr" refers to Dirac, flavor and color matrices. 
An alternative effective potential for composite operator which is a variation of Eq.(|l]) was 



developed in Ref. [jlB] . Here we use the one presented above. 

In the Landau gauge and the one-loop approximation, we can expand the quark two-point 
function as follows 

l G-\ P )=p'-(z s (p 2 ) + t l5 i; P (p 2 )), (2) 

where £ s (p) and £ p (p) are the scalar part and the pseudo-scalar part of the quark dynamical 
mass function, respectively. Substituting (Q) into ([[]) and using the running coupling in the 
two- loop level potential we obtain the following effective potential V — V1+V2 which satisfies 
the renormalization group equation after Wick rotation, 



V l = -2N f N c J 



d*p ^T? s {p 2 ) + T? P {p 2 )+p 2 



(2vr) 4 p 



v 2 = -6N fNc c J [ -pL pEA. 

ic2 JJ (2vr) 4 (2vr) 4 (p - q) 2 
x Mp 2 )Zs(Q 2 ) + Sp(p 2 )E P ( 9 2 ) 

where Nf is the number of the light flavors. This expression is our starting point. The order 
parameter can be determined by the condition 



These extremum conditions for V with respect to £s(p 2 ) and T,p(p 2 ) are equivalent to the 
Schwinger-Dyson equations for the quark propagator in the improved ladder approximation. 
We take £p(p 2 ) = without loss of generality because of the chiral symmetry. We use 
Higashijima-Miransky approximation [5, 6] for the running coupling with an infrared modifi- 
cation |7] 

g 2 (p, q) -> e 2 (max(p 2 , g 2 )), (5) 

where 
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Here, a = 24/ (11 N c — 2nf), rif is the number of flavors which controls the running coupling, 
and pr is a parameter to regularize the divergence at p = 1(Aqcd)- The modified running 
coupling has the same asymptotic behavior as the QCD running coupling, while it gradually 
approaches a constant as p 2 decreases. The value of the running coupling at p = is 2n 2 a/tji, 
where tn = ln(p^) is the infrared regularization parameter. 

As e 2 (max(j9 2 , q 2 )) has no angle dependence in four- dimensional Euclidean space, Eqs.([|) 
and M) are reduced to 



Vi\Z 8 ]=*NfN c j 



d 4 p £|(p 2 ^ 



K 2 [S 5 ] = -12N f N c C 2 J 



(2tt) 4 S|(p 2 )+p 2 
d*p_ E 2 (p 2 ) +p 2 
(2vr) 4 p 2 

d A p e 2 (p 2 ) Y, s {p 2 



-2^/^ln^^, (7) 



(2tt) 4 p 2 E 2 {p 2 )+p 2 
(2tt)* S 2 (g 2 )+g 2 



In this paper we consider the case where the chiral symmetry is exact, i.e. renormalized 
current quark masses are zero. We use the variational method to investigate the vacuum of 
the QCD-like theory and use the following form for Sg(p) as a trial mass function [^]: 

2tt 2 o 



W) = -(qq) 3{p2 + p2 ^ (Hp 2 + p 2 R )r^f , (9) 

where (qq) is the renormalization group invariant vacuum expectation value of qq. This form 
was chosen so that its asymptotic form 

Z S (P 2 ) - -(qq) ^0(\np 2 )^7-\ 

coincides with the one obtained by using the operator product expansion and the renormal- 
ization group equation [14]. This form of dynamical mass function does not break the chiral 
symmetry explicitly, i.e. the axial vector current is conserved^. 



III. EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 



We apply the imaginary time formalism JT7J to calculate the effective potential at finite 
temperature. Since fermion fields obey an anti-periodic boundary condition, we have 



1 We note, however, it was found 15] that the chiral Ward-Takahashi identity is violated unless 
gluon momentum squared is used as the argument for the running coupling. The solution to recover 
the chiral Ward-Takahashi identity in the framework of the improved ladder approximation was 
given in Ref. [16] 
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p A -> o; n = (2n + 1)ttT, (10) 

and the four-momentum integration is replaced with the summation over the Matsubara- 
frequency uo n and the three-momentum integration 

d ' P T± /* (U) 



(2t) 4 (2t) 
in Eqs.(|7|) and Now we have 



Vi = 4N f N c T ]T / 



ri=— oo 



(2tt)3 s|K, P ) + g2 



(2vr) 3 V ^ 

dp e 2 (u; n ,p) £ s (u; n ,p) 



K 2 - -12 N f N c C 2 Tj:j ^ s|K)p) + g , 

xt!4 J; u ::^- d3) 



"(2tt) 3 £ 2 (^ m ,g)+Q 
where Q 2 n = u% + p 2 , Q 2 n = u 2 n + q 2 , and 

^ dq 

4-"»(2tt) 3 

denotes the summation over m and integration over q with the condition Q 2 n > Q 2 m and we 
have replaced e 2 (p 2 ) and £s(p 2 ) with e 2 (u; n ,p) and E^ti^, p), respectively. Their functional 
forms are given below. 

Let us now consider the choice for the running coupling and the quark dynamical mass 
function at finite temperature. We expect that the running coupling decreases as the tem- 
perature increases. Usually a coupling of the form like 

2 (T) 27r 2 q 
9 1 } ln(T 2 /M 2 ) 

is used. However we use a running form 

e 2 (u n ,p) = 2?F " (14) 

ln(^ 2 + p 2 + p 2 R ) 

natural extension of the T = case(Eq.(^)). 
As concerns the dynamical mass function, we use the following function 

S 5 K,p) = -m Kul l" 2 p a 2+p2R) (inU + p' + pl))^" , (15) 

by replacing p 2 with a; 2 + p 2 also in Eq.([5|). Here, (qq) is the variational parameter of the 
effective potential and its value is determined from the minimum of the potential. Note that 
the chiral symmetry is not explicitly broken by the choice Eq.(O). 
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Using the effective potential in Eqs.(|T^) and ( |T3"D with the ansatz Eq.(|14D, Eq.(|15D and 
changing the variables 

t = \n(uj 2 n + p 2 +p 2 R ), 
t n = ln(w 2 +p 2 R ), 



we obtain 



V(a,T) = V 1 (a,T) + V 2 (a,T), 



Vx(a,T) 



AN f N c Ta 2 ~ /■<*> 



(2vr) 2 



r dt Vt 

J tn 



-2t 



t rjtr. 



■IXrX.-T g 



(27T) 
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J t n 
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V2(a,T) 



8X f X c aT 2 a 2 00 ' ^ 



(2^) 2 J=L 
^2 du Ve u — e* 



t -a/2 e -2t 



1 - e*«-* t 2 ~ a (l - e**-*) + a 2 e- 3t 
u i-a/2 e -u 



m ""Tn 



u 2-a(\ _ gtfl-u) _|_ cr2 e - 



3m ' 



(17) 



where a = 2ir 2 a(qq) /3 is the rescaled order parameter. The sum over m, 2 m ', is taken to 
an extent that t > u. 

In order to determine the critical temperature T c , we expand the effective potential as a 
power series of the order parameter with finite coefficients a 2n (T) 



V(a,T) = a 2 (T)a 2 + 0(a A ) 
The coefficient a 2 (T) is calculated as follows 

f rT1 \ dV 
a 2 {T) 



cr 2 =0 



da 2 



2X f X c T 



(2tt) 2 

-4Ta I dt V' 



dt \le l — e tn — — - 

±?„Jtn t 2 ~ a (l-e* 



_n=— 00 

00 



, , > I du v e" — e' m J 

l-e*«"*) 2 tT l-e*«"" 



a/2— lp— u 



The critical temperature T c is determined by the condition, 02 (T c ) = 0. 



(18) 



(19) 



IV. NUMERICAL RESULTS 

A. The Behavior of the Effective Potential at Finite Temperature 

In our calculation we put Nf = 2, X c = 3 and rif = 3. We use the three flavor(nj = 3) 
running coupling since three flavors(u,d,s) contribute to the running coupling in the energy 
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region of the chiral phase transition. Since it has been shown that the quantities such as 
(qq) and f n are not sensitive to the infrared regularization parameter t R at zero temperature 
[jjj, we expect that these quantities are insensitive to t R also at finite temperature as long 
as t R is not too large. In Fig.l we plot the t R dependence of the value of (qq), the location 
of minimum, the pion decay constant f n which is approximately given by the Pagels-Stoker 
formula 



f*=4N r dp Esip) fs (j) E ^m ) 

JlT c 7 (2vr) 4 (S|(p 2 )+P 2 ) 2 V 2 dp 2 J 

and the constituent quark mass m defined by the relation, m = £s(p 2 — 4m 2 ). The fact 
that these values are insensitive to the infrared cutoff at zero temperature has been also 



pointed out by Kugo et al |16] who used the Schwinger-Dyson equation in the improved 
ladder approximation. Therefore we fix the infrared regularization parameter t R to 0.1. 
In Fig. 2 we show that at zero temperature the approximate form of the dynamical mass 
function Sg(p) at the minimum of V well reproduces the exact solution of the Schwinger- 
Dyson equation. Therefore we adopt Eq.(|l5|) as a trial function and use a as a variational 
parameter at finite temperature. We have confirmed that the quantities such as (qq)r and 
f-w(T) are insensitive to the infrared regularization parameter tR at finite temperature. 

The effective potential at finite temperature is shown in Fig. 3. At zero temperature 
the chiral symmetry is spontaneously broken. As temperature grows the minimum of the 
potential continuously approaches zero. Hence we have second order phase transition at 
T = 0.185(Aq(7£)). To determine the value of Aq C d we use the Pagels-Stoker formula and 
fix the pion decay constant / w to 93 MeV at zero temperature and have Aqcd = 737MeV. 
Then the restoration of the chiral symmetry occurs at T c = 136MeV. 

The temperature dependence of the order parameter (qq)r is shown in Fig.4. As concerns 
the temperature dependence of (qq)r we refer, for example, to that of the free pion gas 
approximation: 

<w>o s/r ( } 

This behavior is also shown in Fig.4 and there is a slight difference between our result and 
Eq.(P0D in the low temperature region. 

Once the quark dynamical mass function is determined from the minimum of the effective 
potential, we can calculate the pion decay constant /„. at finite temperature using the Pagels- 
Stoker formula extended to finite temperature: 

I IT\ 2 - iN T V / dp S sK,p) (y ,, , „\ Qi <EsK,p) \ 

u( t) - ^ ( ^ KiP) + <K) » (^-p) - - -^r) ■ 

where Q\ = cj 2 + p 2 . The behavior of /^(T) is shown in Fig.5. 



B. Critical Exponents 

Since we have a second order phase transition, we can estimate the critical exponents. 
To estimate the critical exponents we define the reduced temperature r and exponents j3, 7 
and (3' as follows: 
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\(QQ)t\ 


' ? 


« 2 (T)| 









(T < T c ), 



(21) 

In the spin system the exponent 7 is defined for the magnetic susceptibility x{T)- In 
case of the chiral phase transition, x{T) is defined as 



X(T) - m)T 



dmn 



mo=0 



where mo is the current quark mass. Our exponent 7 agrees with that of the magnetic 
susceptibility, if the effect of the finite current mass near T c is represented as 

V(cr, Too, T) ~ a 2 (T)a 2 + b(T)m a, 

and, furthermore, if b(T) does not have a singular behavior. 

In Landau theory of second order phase transition, the assumption that 7 = 1 leads to 
(3 = 1/2. In Ref. [10] they study the critical exponents within the Landau theory. In general 
7 deviates from 1. So we first examine the critical exponent 7 for three cases of the infrared 
regularization parameter tn = 0.05,0.1,0.15. We determine the exponent 7 from a slope of 
the linear log fit 

ln|a 2 (T)| = 7 lnr + C, 

where C is independent of r. In order to extract the exponent we use the \ 2 fitting(see 
Fig. 6) and have 

1.032 ±0.006 
7 = I 1.008 ±0.001 
1.041 ±0.008 

These results are consistent with the Landau theory. 

In the same way(see Fig. 7) we determine the exponents (3 for £r = 0.1 and have 

(3 = 0.455 ±0.008, 

which is appreciably different from the Landau theory. 

Next, let us estimate (3'. We used the Pagels-Stoker formula and Eq.flTBD for calculating 
fn{T), thus we have 

U{T) = \{qq) T \G{{qq) T ,T). 

We have checked numerically that G((qq)T, T) is free from the singular behavior in the limit 
t — >• 0. It follows that the behavior of f n (T) near T c is the same as that of the condensate: 

U(T) ~ (22) 





(0.05\ 


for t R = 


0.10 




\ 0.15 / 
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Thus we have /?' = /? = 0.455. The same analysis can be repeated for the exponent of the 
mass gap T*s{ n — 0, p = 0). As far as we use Eq.(15) as a trial mass function, the exponent 
of the mass gap agrees with the exponent (3. 

Note that we use Eq . (plf) to evaluate the critical exponents. As a rule, for D > 4 (D = 
dimension of space-time) one obtains the universal prediction of the mean field theory. How- 
ever, for D < 4 the infrared divergences which appears in the coefficients of the power series 
expansion invalidate the predictions of the mean field theory []I9" |. Here we are interested in 
the case of D = 4, which is critical. Although our results are well fitted by Eq.fl2l~D, there 
may be logarithmic corrections to the power behavior assumed in Eq.( 



V. SUMMARY AND DISCUSSION 

In summary, we studied the chiral phase transition at finite temperature in the QCD-like 
gauge theory. We made use of the CJT effective potential in the improved ladder approxi- 
mation and took account of the logarithmic behavior in the quark dynamical mass function. 
In QCD-like gauge theory it has been pointed out that at zero temperature the physical 
quantities such as (qq), fx and m are quite stable under a change of the infrared cutoff. We 
confirmed that the temperature dependence of these quantities are also stable under a change 



of t R (see also Ref. |TT|). Then we put t R = 0.1 and investigated the temperature dependence 
of the effective potential and found the second order phase transition at T c = 0.185Aqct>. 
We fixed the value of the pion decay constant to 93MeV at zero temperature and had 
Aqcd = 737MeV, therefore we obtained T c = 136MeV. We also examined the temperature 
dependence of (qq)x and f n (T). 

As we have the second order phase transition, we calculated the critical exponents 7 
and (3 which are the exponents of the quadratic term of the effective potential and the 
order parameter, respectively. We examined 7 for three cases of the infrared regularization 
parameter and found 7 ~ 1. This result is consistent with the Landau theory. However we 
found that f3{= 0.455 ± 0.008) appreciably deviates from the value of the Landau theory 
1/2. Our approach is somewhat related to the mean field approximation. However we have 
introduced the logarithmic behavior in the quark dynamical mass function. This may be one 
of the main reasons that (3 deviates from the Landau theory. Our results are different from 
0(4) Heisenberg model. We also calculated the exponents of f n (T) and T>s( n = 0,p = 0) 
and they coincided with (3. 

In calculating the temperature dependence of several quantities related to the chiral phase 
transition and estimating the critical exponents, the effective potential and/or the Schwinger- 



Dyson equation with a variety of approximation have been used[9, 10,11, 12]. In Ref. JTU] they 
have used a variety of CJT effective potential, which is modified by using Schwinger- Dyson 
equation in Eq.(jl]), and the trial mass function without logarithmic behavior. However they 
have discussed the critical phenomena within a Landau theory. In Ref. || they have used 
the Coulomb gauge Schwinger-Dyson equation with an instantaneous approximation and 
have found (3 = 0.49. In recent work based on the confining Schwinger-Dyson equation 



14 , it has been found that (3 = 0.46 ± 0.04. On the other hand, in Ref. |TT| they have 
used the Schwinger-Dyson equation in the improved ladder approximation and have found 
that (3 = 0.171. Although their model is in the same framework as ours, however, this 
value largely deviates from the Landau theory and our result. We calculate the effective 
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potential by using the trial mass function, while the authors of Ref. JTT| solved the Schwinger- 
Dyson equation. Furthermore, there are important differences in the handling of the gluon 
propagator at finite temperature and in the choice of the argument of the running coupling, 
they used g 2 (p, q) = g 2 (p 2 + q 2 )- 

In this paper, we have obtained a large value of Aqcd- This may be partly due to the 
fact that we did not consider the effect of the color confinement. In fact, in Ref. |20|| , they 
have used an observed value of Aqcd to reproduce the pion decay constant. 

As mentioned in Sec. IV, some critical exponents are defined in relation to the current 
quark mass. Therefore, it is very interesting to generalize our model to the case of a finite 
current mass and to study the critical exponents and the scaling relations. 
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Figure Captions 



Fig. l:t R dependence of /„-, |(gg)| 1//3 ,m at T = 0. The mass unit is Aqcd- 
Fig. 2:The dynamical mass function of quarks for t R = 0.1, T = 0. The curves correspond 
to the solution of Schwinger-Dyson equation (solid line) and the trial mass function Eq.(15) 
at the minimum of V (dashed line). 

Fig. 3:Effective potential at finite temperature as a function of Kqq)] 1 ^ 3 for case t R = 0.1, 
Nf = 2 and N c — nf — 3. The curves show the cases T = (solid line), T = 0.15(Aq C d) 
(dashed line), T = 0.25(A QC d) (dotted line). 

Fig. 4:The temperature dependence \{qq)r\ for t R = 0.1. The dashed line corresponds to 
the free pion gas approximation. 

Fig. 5:The temperature dependence of fn{T) for t R = 0.1. 

Fig. 6:Linear log fit to In |ei2(T)| for £r = 0.05(triangle), O.l(cross), 0.15(circle) respectively. 
Here the absolute normalization of a 2 (T) is taken arbitrary. 
Fig. 7:Linear log fit to \na min for t R = 0.1. 
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